Condensation dynamics in a quantum-quenched Bose gas 
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By quenching the strength of interactions in a partially condensed Bose gas we create a "super- saturated" 
vapor which has more thermal atoms than it can contain in equilibrium. Subsequently, the number of condensed 
atoms (A^o) grows even though the temperature (T) rises and the total atom number decays. We show that 
the non-equilibrium evolution of the system is isoenergetic and for small initial A^o observe a clear separation 
between T and A'o dynamics, thus explicitly demonstrating the theoretically expected "two-step" picture of con- 
densate growth. For increasing initial A'o values we observe a crossover to classical relaxation dynamics. The 
size of the observed quench-induced effects can be explained using a simple equation of state for an interacting 
harmonically-trapped atomic gas. 
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Non-equilibrium dynamics of interacting quantum systems 
are generally far less understood than the corresponding equi- 
librium many-body states 1 1 1. Of particular interest are many- 
body dynamics of both the order parameter and the excita- 
tions in a system close to a phase transition. From a the- 
oretical point of view, a clean and well defined way to in- 
duce and study non-equilibrium quantum dynamics is a rapid 
"quantum quench" | 2 1 of a single Hamiltonian parameter. Ul- 
tracold atomic gases are very well suited for such quantum 
quench experiments. In addition to the possibility to dynam- 
ically vary microscopic Hamiltonian parameters, they feature 
near-perfect isolation from the environment and characteris- 
tic many-body timescales (ranging from milliseconds to sec- 
onds) that are experimentally resolvable and allow real-time 
non-equilibrium studies. 

In this Letter, we introduce a quantum quench of the inter- 
action strength in an atomic Bose gas as a tool to study the 
dynamics of Bose-Einstein condensation ISHTSll. Earlier ex- 
periments highlighted the importance of bosonic stimulation 
in condensate formation ifTOl , but could not quantitatively ad- 
dress the theoretically debated interplay of energy redistribu- 
tion and coherence development in the system (SHSl . The use 
of a quantum quench of the interaction strength allows us to 
study these two processes in parallel. The quench induces 
a growth of the condensed atom number in a degenerate gas 
without any removal of thermal energy; we explain this effect 
with a simple theoretical model and experimentally study its 
real-time dynamics. We explicitly show that the post-quench 
non-equilibrium evolution of the system is isoenergetic, and 
directly reveal the theoretically postulated "two-step" picture 
of condensation |4-7J . As expected, close to the critical point 
the growth of the condensed atom number lags behind the en- 
ergy redistribution in the thermal component of the gas. Mov- 
ing away from the critical point, we also observe a crossover 
to effectively one-step condensation dynamics governed by a 
classical relaxation process. 

In an ideal Bose gas the number of condensed atoms, A^o. 
depends only on the total atom number A^ and the temperature 
T. In a partially condensed cloud at a given T the number of 
atoms in the thermal component, N\ is saturated at the critical 



value for condensation, Nc{T), and we have A^o = ^ — ^c- 
In experiments on harmonically trapped atomic gases the in- 
teractions, characterized by the s-wave scattering length a, 
change this picture in two ways (see Fig. [T]). First, they in- 
duce a shift of the critical value Nc{a^T), which was accu- 
rately measured in |16|. Second, the thermal component is 
not saturated - the presence of the condensate allows A^' to 
grow above Nc |17|. Taking these effects into account, near 
the critical point we can write the equation of state for an in- 
teracting atomic gas in thermal equilibrium 1 18 J : 

N = N,^ SoN^^^ + No . (1) 
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Here SqNq^ is the additional number of atoms accommo- 
dated in the thermal component due to non- saturation ef- 
fects, i.e., A^' = A'c + SqNq^^. The non-saturation coeffi- 
cient So (X a'^^^T'^ can be calculated using mean-field theory 
El [Hi. (In our experiments 5'o 10^.) 

As illustrated in Fig. [T] the fact that the equilibrium num- 
ber of condensed atoms (at a fixed A^ and T) depends on the 
strength of interactions opens the possibility to use a quantum 
quench of the scattering length to induce non-equilibrium A'o 




FIG. 1: (Color online) Inducing non-equilibrium dynamics by a 
quantum quench of the interaction strength. The equilibrium con- 
densed atom number A'o, calculated according to Eq. ([T]), is plotted 
versus the total atom number A^ for a fixed temperature T — 200 nK, 
our trapping parameters, and different scattering lengths a. The red 
arrow indicates the direction of the quench. The three absorption 
images of atomic clouds released from the trap show the growth 
of the condensate (appearing in red) over 1 s following an 
a = 275 62 ao quench. 



2 



dynamics. Here we plot the solutions of Eq. ([T]) for a fixed T 
and different scattering lengths a, where ao is the Bohr radius. 
If the system is prepared in equilibrium at a high scattering 
length and then a is quenched to a lower value the gas be- 
comes super-saturated, having too large N' and too small Nq. 
Consequently Nq must grow even without any active cooling 
of the gas. Compared to the classical quench of the thermal 
energy 1 10-14 L the quantum quench of a has the advantage 
that it does not directly affect N and the initial thermal oc- 
cupations of the single-particle excited states. This allows us 
to study in parallel the induced non-equilibrium evolution of 
both the coherence and the quasi-thermal energy distribution 
in the system. 

For our experiments we use an optically trapped cloud of 
^^K atoms in the |F, m^) = |1,1) hyperfine ground state 
|[T9l , in which a can be tuned via a Feshbach resonance cen- 
tred at 402.5 G f20l. The geometric mean of the harmonic 
trapping frequencies in our nearly isotropic trap is cD/27r ^ 
75 Hz, the temperature of our clouds is T ^ 200 nK, and the 
total atom number TV ^ (3 - 4) x 10^. 

We initially prepare an equilibrium cloud just below the 
condensation temperature at a scattering length a* = 275 ao 
and then rapidly quench the scattering length to a lower final 
value, , by changing the externally applied magnetic field 
over Tq = 10ms. After the quench we follow the evolution 
of the system for times up to t = 10 s, extracting A'o. N and 
T from absorption images taken after 1 8 ms of time-of-flight 
(TOF) expansion from the trap | 21 1. An example of such a 
data series, with = 62 ao and an initial condensate atom 
number A^q ~ 4 x 10^, is shown in Fig. [2] We collected 
a total of 20 such experimental series, each with a different 
combination of in the range 52 — 97 ao and A^q in the range 
(1-40) X 10^. 

The quench time Tq and the values simultaneously sat- 
isfy several requirements: 

(i) Even for fixed A'o the spatial size of the condensate de- 
pends on a. Our Tq is sufficiently long compared to the trap 
time (cjTq ^ 5) to allow adiabatic adjustment of the conden- 
sate shape at the initial A^o |22|. 

(ii) At the same time the quench must be diabatic with re- 
spect to the exchange of particles between the thermal cloud 
and the condensate, so we choose values small enough for 
the elastic collision rate 7ei 1 23 1 to be much smaller than 1/rq. 

(iii) The values are large enough for the system to 
converge towards new equilibrium at long times, rather than 
forever remaining in an intrinsically non-equilibrium state 
MM- 

As shown in Fig. |2] the number of condensed atoms A'o 
clearly grows following the quench, before eventually decay- 
ing at much longer times due to the mundane reasons of slow 
A^ decay and a background heating rate of about 1 nK/s 1241 . 
Since at short times A'o grows while A' decays and T rises, 
this increase in the number of condensed atoms is unambigu- 
ously an interaction effect. Interestingly, the temperature also 
shows a fast initial rise which is clearly associated with the 
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FIG. 2: (Color online) Non-equilibrium dynamics following an in- 
teraction quench from a* = 275 ao to = 62 ao. The condensed 
atom number Ao, total atom number A, and temperature T are plot- 
ted versus the time t after the quench. Each data point is an average 
of 3 - 7 experimental shots and the error bars are statistical. Ao first 
grows towards the new equilibrium and then slowly decays due to the 
slow A decay and a small background heating rate of 1 nK/s. The 
A decay has a characteristic timescale of 100 s and is essentially 
linear over 10 s. In addition to the constant background heating rate, 
the temperature shows a fast initial rise induced by the quench. The 
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inset shows the evolution of Aq ^ , which decays linearly at long t. 
The dashed red lines show fits to the long-t data, used to extract Aq 
and (see text). 



interaction quench. 

In order to quantitatively study the non-equilibrium effects 
that occur on short timescales (^ 1 s) after the quench, we 
need to eliminate from our analysis the long-term (quasi- 
static) drifts of A" and T. Specifically, we need to experi- 
mentally extract the "target" final values A^q and that the 
system would tend to in absence of the slow background heat- 
ing and atom number decay (see dashed red lines in Fig. [2]). 
The target temperature is simply determined by subtract- 
ing the constant background heating rate, but the extraction of 
A'q is a bit more subtle. The slow linear decay of A" and rise 
of T both to leading order correspond to a linear decrease of 

2 /5 

N — Nc with time. Further, for small condensates, the SqNq^ 
term in Eq. ([T]) is significantly larger than A'o. Hence, at long 
times A^q'^^ decays linearly, allowing us to extract (A^q )^/^ by 
linear extrapolation to t = 0. 

We qualitatively anticipated the quench-induced growth of 
A'o by solving Eq. ([T]) under the constraint of fixed T (see 
Fig. [T]). However we can also understand why the quench 
must lead to an increase in the temperature of the cloud. For 
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FIG. 3: (Color online) Isoenergetic non-equilibrium evolution, (a) 
For a quench from = 200 nK, and a* = 275 ao to a-^ — 62 ao, 
the three lines show the calculated A^o (bottom, black), Nq in the 
isoenergetic picture (middle, red) and in the isothermal picture 
(top, blue), (b) Measured AA^o, for 20 experimental series with var- 
ious A^o and values, is plotted versus the predictions of the isoen- 
ergetic model. The liner fit to the data (solid blue line) gives a slope 
of 0.9 ± 0.1 and an offset of (3 ± 1) x 10^ 



our experimental parameters the average energy of thermal 
atoms is approximately S/c^T ^ h x 13 kHz, while the en- 
ergy of the condensed atoms (including kinetic, potential, and 
interaction energy at a^) is less than h x IkHz. Therefore, 
during the non-equilibrium evolution of the system after the 
quench, the atoms moving from the thermal cloud into the 
condensate take with them much less than "their share" of en- 
ergy. Hence an isothermal non-equilibrium evolution would 
not conserve energy. To explicitly test whether the evolution 
is isoenergetic, we make a simple prediction that relates the 



A^Q under the 



changes AT = - T' and AA^o = N/ 
constraint of constant energy. In analogy with the standard 
models of evaporative cooling [ 25 , 26], we get that to leading 
order (for small condensed fractions) the small fractional in- 
crease in T should be equal to the increase in the condensed 
fraction: 



AT 



ANq 
N 



(2) 



For our 20 data series with different Nq and values we 
observe AT/T' = (1.2 ± 0.3)ANo/N, in good agreement 
with the prediction of Eq. ([2]). We can now also predict ANq 
for any Nq and , by numerically solving Eq. ^ under the 
constraint set by Eq. ^ (see Fig. [3]). As shown in Fig.[3jb), the 
measured AA^o follows our theoretical predictions, although 
we systematically observe slightly smaller Nq increase than 
predicted |27|. 

In Fig.|4]we compare and contrast the evolution of T and 
Nq during the system's approach to the new equilibrium, for 
the same 20 data series shown in Fig.[3jb). 

The temperature exhibits classical relaxation dynamics, i.e., 
T exponentially approaches on a timescale tt that de- 
pends only on 7ei. As shown in Fig.Qa), we find that tt cor- 
responds to 2.6 collisions per particle (see also, e.g., 12911301 ). 
We observe no dependence of tt on A^q. 

The A^o dynamics are more intriguing. For Nq to grow, two 
conceptually distinct steps must take place: (1) the redistribu- 
tion of the kinetic energy within the gas (seen in the T dynam- 




FIG. 4: (Color online) Condensation dynamics, (a) The T dynamics 
agree with the classical relaxation picture; the linear fit (red line) 
gives a relaxation time tt = (2.6 di 0.1)/7ei. (b) In the bosonic 
stimulation picture Ro oc Fb (see text). The linear fit (solid blue 



line) to \ow-Ro data gives Rq 



(B) 



(9 ± l)rB. (c) One- versus 



two-step condensation. Solid red line shows the one-step classical 
relaxation prediction, Rq . The blue shaded area, corresponding to 
the blue line in (b), shows the two-step bosonic stimulation result, 
R^^\ For increasing Rq we observe a crossover between the two 
types of dynamics. The inset shows the same data plotted on a linear 
scale. 



ics) and (2) the merging of the accumulated low-energy atoms 
into the coherent condensate. In this two-step picture, (only) 
the second step depends on A^q due to bosonic stimulation, 
which enhances scattering into an already highly occupied 
state. Specifically, the initial A^o growth rate, Rq = No{t = 
0), should be proportional to Fb = A^q 7ei ^m/(^b^), where 
A/1 ^ c?!^ A{Nq^^) is the difference between the initial and 
final chemical potentials I 8H10I . 

However, the two-step picture is experimentally relevant 
only if the second step is the slower, "rate limiting process", 
so that A^o growth lags behind the energy redistribution. If the 
first step is the rate limiting process, and the particles scatter- 
ing to low energies essentially immediately join the conden- 
sate, then we effectively have a one-step process, described 
by the classical relaxation model. In that case we expect Rq 
to be simply given by R^l"^^ = ANq/tt ^ ANq 7ei/2.6. We 
qualitatively expect a crossover from the bosonic stimulation 
(two-step) to the classical relaxation (one-step) behaviour as 
we move away from the critical point by increasing A^q. In 
essence, Rq should always reflect the slower of the two pro- 
cesses. 

In Fig. [4jb) we plot Rq versus Fb EH. The low Rq data 
shows the expected proportionality and the linear fit (solid 
blue line) gives Rq ^ = (9±1)Fb. For large Fb values we see 
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a systematic downwards deviation from this fit. This is what 
is expected once Rq ^ exceeds the rate at which the classical 
relaxation of the thermal component feeds atoms into the low- 
energy states. In this regime the experimentally observed Rq 
should be lower than R^"^ and given by i?Q^^^ . 

In Fig. |4jc) we directly compare the two-step and one- step 
pictures, and show the crossover between the two types of be- 
haviour. Here we re-plot the Rq data versus AA^o 7ei, so that 
the one-step classical result, Rq , corresponds to the straight 
line shown in red. The bosonic stimulation result, Rq % is 
now not a universal curve, since it depends on a different set 
of parameters; the blue line in Fig. Qb) here maps onto the 
blue shaded area. We now explicitly see that the small Rq 
data lies systematically below the classical relaxation predic- 
tion, as expected in the two-step picture. However we also see 
that the data is consistent with both theories in the crossover 
region where R^"^ ^ i?Q^^^ , and eventually agrees better with 
Rq ^ for the largest Rq values. 

In conclusion, we have used a quantum quench of the in- 
teraction strength to create a super-saturated non-equilibrium 
Bose gas and study its dynamics. We have shown that the 
non-equilibrium evolution of the system is isoenergetic and 
that the quench-induced changes in the condensed fraction 
and temperature of the gas can be accounted for using a sim- 
ple equation of state for an interacting gas in thermal equilib- 
rium. Moreover, for the first time we directly compare and 
contrast the energy-distribution and coherence dynamics, and 
clearly resolve the two theoretically expected steps in the con- 
densation process. Here we focused on the case of small but 
non-zero initial condensates; with this case understood, in the 
future it should be possible to use a similar quantum quench 
to drive the system through the critical point. In that case, in 
absence of the initial condensate "seed", it should be possible 
to study the stochastic effects associated with the spontaneous 
symmetry breaking and initial condensate formation. 

We thank R. Fletcher and A. Gaunt for comments on the 
manuscript and N. Tammuz for experimental assistance. This 
work was supported by EPSRC (Grants No. EP/G026823/1 
and No. EP/IO 105 80/1) and a grant from ARO with funding 
from the DARPA OLE program. 



[1] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, 

Rev. Mod. Phys. 83, 863 (2011). 
[2] R Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 (2006). 
[3] D. W. Snoke and J. P Wolfe, Phys. Rev. B. 39, 4030 (1989). 
[4] H. T. C. Stoof, Phys. Rev. Lett. 66, 3148 (1991). 
[5] Y. Kagan, B. V. Svistunov, and G. V. Shlyapnikov, Sov. Phys. 

JETP 75,387 (1992). 
[6] D. V. Semikoz and 1. 1. Tkachev, Phys. Rev. Lett. 74, 3093 

(1995). 

[7] See Chapters by Yu. Kagan and H. T. C. Stoof in Bose- 
Einstein Condensation, edited by A. Griffin, D. W. Snoke, and 



S. Stringari (Cambridge University Press, Cambridge, 1995). 
[8] C. W. Gardiner, P Zoller, R. J. Ballagh, and M. J. Davis, Phys. 

Rev. Lett. 79, 1793 (1997). 
[9] C. W. Gardiner et al, Phys. Rev. Lett. 81, 5266 (1998). 

[10] H.-J. Miesner et al. Science 279, 1005 (1998). 

[1 1] M. Kohl et aL, Phys. Rev. Lett. 88, 080402 (2002). 

[12] S. Ritter et al, Phys. Rev. Lett. 98, 090402 (2007). 

[13] M. Hugbart et aL, Phys. Rev. A 75, 01 1602 (2007). 

[14] C. N. Weiler et al. Nature (London) 455, 948 (2008). 

[15] M. C. Garrett et aL, Phys. Rev. A 83, 013630 (2011). 

[16] R. P Smith, R. L. D. Campbell, N. Tammuz, and Z. Hadzibabic, 
Phys. Rev. Lett. 106, 250403 (2011). 

[17] N. Tammuz etal, Phys. Rev. Lett. 106, 230401 (2011). 

[18] R. P Smith and Z. Hadzibabic, arXiv: 1203.2063 (2012). 

[19] R. L. D. Campbell etaL, Phys. Rev. A 82, 063611 (2010). 

[20] M. Zaccanti et aL, Nature Physics 5, 586 (2009). 

[21] In the initial non-equilibrium regime, at short times after the 
quench, T is strictly speaking not defined. However, in our ex- 
periments the quench results in changes of the apparent tem- 
perature by < 5%. For such small deviations, at all times an 
equilibrium thermal distribution fits the data very well and pro- 
vides a good measure of the thermal energy. For fitting the ef- 
fective temperature we exclude the central portion of the cloud 
corresponding to one thermal radius. 

[22] M. R. Matthews etaL, Phys. Rev. Lett. 81, 243 (1998). 

[23] 7ei = V2{n)vth(yo, where (n) is the average density of thermal 
atoms at the critical point, ^;th = i/S/cBT/Trm and ctq = Stto^. 
We calculate (n) using an ideal Bose gas distribution at the crit- 
ical point. 

[24] The main sources of the slow atom loss and the background 
heating are off-resonant spontaneous scattering of photons from 
the trapping laser beams and small fluctuations in the trapping 
laser power. The one -body loss due to collisions with the back- 
ground gas in the vacuum chamber has an even longer charac- 
teristic timescale of ^ 400 s. 

[25] O. J. Luiten, M. W. Reynolds, and J. T. M. Walraven, Phys. Rev. 
A 53, 381 (1996). 

[26] W. Ketterle and N. J. van Druten, in Advances in Atomic, 
Molecular, and Optical Physics, edited by B. Bederson and H. 
Walther (Academic Press, San Diego, 1996), Vol. 37, pp. 181- 
236. 

[27] Since our simple model does not include any non-linear decay 
process, this small discrepancy could be due to a small three- 
body loss of condensed atoms at short times after the quench. 
While we do not observe any three-body effects in a thermal 
gas, A^o could be slightly affected due to the higher density of 
the condensate. For A^o ^ 10^ the loss rate per atom is more 
than an order of magnitude larger for the condensed than for 
the thermal component Il28|. Another possible reason for the 
discrepancy is that the non-zero quench time Tq results in an 
effective a* that is slightly lower than 275 ao- 

[28] J. Soding et aL, Appl. Phys. B 69, 257 (1999). 

[29] C. R. Monroe et aL, Phys. Rev. Lett. 70, 414 (1993). 

[30] B. DeMarco et aL, Phys. Rev. Lett. 82, 4208 (1999). 

[31] To experimentally extract Rq we heuristically use an exponen- 
tial fit to the short-t data, taking into account the slow long-term 
decay of A^o- In the bosonic stimulation picture this functional 
form is not exact, but within experimental precision it fits the 
data with non-zero A^o sufficiently well. For example, the solid 
blue line in Fig.[2]is based on such a fit. 



